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ABSTRACT 


Various Symbolic Manipulation Programs have been tested to check the 
functioning of their commands and suitability under various operating 
systems. Support services for SMP have been found to be relatively better 
than the one for MACSYMA. The graphics facilities for MACSYMA does not work 
as expected under UNIX operating system. Not all the commands for MACSYMA 
function as described in manuals. 

Shape representation has been a central issue in computer graphics and 
computer-aided design. Aside from appearance, there are other application 
dependent, desirable properties like continuity to certain order, symmetry, 
axis- independence and variation-diminishing properites. Several shape 
representations are studied, which include the Osculatory Method, a 
Piecewise Cubic Polynomial Method using two different slope estimates. 
Piecewise Cubic Hermite Form, a method by Harry McLaughlin, and a Piecewise 
Bezier Method. They are applied to collected physical and chemical data. 
Relative merits and demerits of these methods and examined. 

Kinematics of a single link, non-dissipative robot arm is studied using 
MACSYMA. Lagranian is set-up and Lagrange's equations are derived. From 
there, Hamiltonian equations of motion are obtained. Equations suggest that 
bifurcation of solutions can occur, depending upon the value of a single 
parameter. Using the characteristic function W, Hamilton- Jacobi equation is 
derived. It is shown that H-J equation can be solved in closed form. 
Analytical solutions to H-J equation are obtained. 



I. INTRODUCTION 


This project started in August, 1983 and ended in December, 1986. From 
August, 1983 till August, 1986; Dr. Shantilal Shah was the Principal 
Investigator. Dr. Shah left Hampton University in September, 1986. Dr. 
Alkesh Punjabi was Principal Investigator from September, 1986 to December, 
1986. Dr. John Wiggs, Dr. Maria Lam and Dr. Alkesh Punjabi have been 
Associate Investigators in this project for varying durations. Dr. John 
Wiggs is no more with Hampton University. 

Hampton University has benefited by this project in a number of ways. 

A number of students learnt how to use symbolic manipulation programs like 
MACSYMA, SMP, NUMATH and others. These students were black minority 
students. The working knowledge and experience of symbolic manipulation 
programs and their application to problems will definitely help these 
students when they enter the market place for job. 

The faculty members who were involved in the project have also 
benefited. They were able to learn and use one of the most important tools 
in modern mathematical and physical sciences research. It is generally 
accepted that symbolic manipulation programs are such a tool. This project 
helped the faculty in carrying out research using MACSYMA and SMP. 

This project has also been helpful to the recently started graduate 
program in Applied Mathematics by way of the graduate research 
assistantships made available to graduate students in this program. The 
work done under this project may be the dissertation material for one 
candidate for M.S. in Applied Mathematics. The computer, the hardware and 
software made available to Hampton University under this project has been of 
immense value to the University and especially to the Computer Science 



Department. These computational capabilities have facilitated the teaching 
of a number of courses in the computer science curricula. It has also 
helped the research activities of mathematics, computer science, and physics 
departments. 

In short, one can definitely assert that this project has been of great 
benefit and help to the University in many ways. 

On the other hand, NASA LaRC has also benefited by this project. NASA 
LaRC was provided access to MACSYMA and SMP on VAX computers at Hampton 
University. Also, under this project, the MACSYMA, SMP, and other programs 
were tested and the defected commands and functions which did not work or 
did not work as postulated were detected and pointed out. It is hoped that 
these aspects of the project have been of use to NASA LaRC. 

Three contributed papers were presented on research done under this 
project. The details of these are given elsewhere in this report. 

The rest of this report is organized as follows: 

Sections 2, 3, and 4 detail the work or research done by Dr. Shah, Dr. 
Lam and Dr. Punjabi respectively under this project. Each section is 
written by the faculty member concerned. Section 5 gives the details of the 
contributed papers. 
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II. Evaluation of Some Symbolic Manipulation Software 


In 1984, Hampton University purchased VAX-11/780 computer with 4 
megabytes memory with VMS operating system and installed REX-MACSYMA 
software purchased from Symbolics Inc. MACSYMA is a symbolic manipulation 
software which was developed at M.I.T. over the period of the last 20 years. 
It was developed and resided on Digital Equipment Corporation PDP-10 
Computer. Many scientists at M.I.T. had also created '1 ibr ary packages to 
solve various problems in pure and applied mathematics and it was made 
available for public use under shared libraries under MACSYMA. MACSYMA has 
also capabilities for point-graphs and 2 and 3 dimensional line-graphs. 
MACSYMA software became an important tool in research and development in 
solving difficult problems in the fields of pure and applied mathematics, 
engineering, physics, aerodynamics and several other disciplines. In 1983, 
M.I.T. sold the marketing and support of MACSYMA to SYMBOLICS Inc. 

SYMBOLICS Inc. started distribution of MACSYMA software which ran on 
VAX-11 /780 computer under VMS operating system in 1983. When it was 
installed on VAX-11/780 computer running under VMS at Hampton University in 
1984, it was able to support most of MACSYMA commands available in the 
language. The investigators spent spring and summer of 1984 in beta testing 
of the software. The line-graphic support was not available on REX-MACSYMA. 
There were some other bugs in few MACSYMA commands. The shared library 
packages were available on REX-MACSYMA software for use, but none of them 
were working. This information was conveyed to Symbolics Inc. at the 
MACSYMA Conference at Schnectudy, N.Y. in summer of 1984. People from 
Symbolics Inc. told the investigators that shared libraries were not the 
part of MACSYMA software package, they were not responsible for its support, 
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and Hampton University was not supposed to have it on its machine. The 
Symbolics Inc. also acknowledged that they had problems in supporting of 
line-graphics, both two and three dimensions, capabilities of MACSYMA at the 
time and they were working on the problems. The execution speed of various 
commands of REX-MACSYMA was also considerably slower compared to the 
software running at M.I.T. 

In spring of 1984, Wollongong Inc. came out with a software named 
EUNICE, which works like UNIX ope'rating system. It supports all of the UNIX 
operating system commands and runs under VMS operating system. When it is 
installed on VAX-11 machine running under VMS operating system, the user has 
advantage of using both VMS and UNIX operating system commands. EUNICE 
software imulates UNIX operating system running under VMS. Investigators 
inquired with Symbolic Inc. whether they had developed a version of MACSYMA 
which runs under EUNICE operating systems and whether EUNICE-MACSYMA will 
support line-graphics capabilities of MACSYMA. Symbolics Inc. informed the 
investigators that they had EUNICE-MACSYMA and the version supported line 
graphics. In 1985, Hampton University investigators purchased and installed 
EUNICE operating system from Wollongong Group Inc., and swapped REX-MACSYMA 
for EUNICE-MACSYMA with Symbolics Inc. Investigators also purchased SMP 
software developed by Inference Corporation. SMP is a symbolic manipulation 
program like MACSYMA, which has capabilities to perform symbolic 
computations and graphics for two and three dimensions. SMP software can 
run under either VMS operating system or UNIX operating system on VAX-11 
computers. Investigators have performed comparison of execution times of 
SMP commands and MACSYMA commands running under VMS operating system. SMP 
software supports graphics capabilities on various graphics terminals. 
EUNICE-MACSYMA, like REX-MACSYMA can do point plots, but does not support 
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line graphics either on Digital's GIGI graphics terminals or Tektronix 
graphics terminals. This was again brought to the attention of people at 
Symbolics Inc. Inquiries were again made with Symbolics Inc. whether the 
UNIX version of MACSYMA available from Symbolics Inc. was able to support 
line-graphics and whether they can provide the names of the institutions or 
companies which are currently users of UNIX-MACSYMA. The reply of Symbolics 
Inc. was that UNIX-MACSYMA does support line graphics on Tektronix terminals 
and due to some Symbolic Inc. policy, they cannot give out the names of 
UNIX-MACSYMA users. Hampton University installed UNIX operating system on 
its VAX-1 1/780 computer in February of 1986 and also installed UNIX version 
of MACSYMA on VAX-1 1/780 computer. Experiments were performed to test 
graphics capabilities of UNIX-MACSYMA on tektronix terminals. The results 
of these experiments showed that, the software was not working correctly as 
desired. Some of the drivers routines required to support graphics were not 
available. Symbolics Inc. tried to evade the responsibilities for not 
having driver routines instead they tried to find fault with UNIX operating 
system and tektronix graphic terminals. Finally in June of 1986, Jeffery P. 
Golden at Symbolics Inc., looked into the problem and found out that 
graphics capabilities on UNIX-MACSYMA were not working properly. In August 
1986, at the Stanford University meeting, he told that he and other people 
at Symbolics Inc. were working on the problem and it will be sometime before 
the problem is corrected. 

Since Symbolics Inc. acquired the distribution of MACSYMA software, 
their support for the software is totally unsatisfactory and the graphic 
capabilities (line) is not working on any of their versions REX, EUNICE or 
UNIX. So far, they are not truthful about giving information on graphics 
capabilities of MACSYMA. They are charging about $350.00 for support 
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service on MACSYMA, and customers do not get any support from the service. 
The shared libraries are not available for MACSYMA users. 

SMP software is comparatively new and young in the market. Some of its 

commands are also not working properly. But their support is very 

excellent. The graphics capabilities on SMP is very good and works 

correctly. Some of the commands like 'XTHRIT which is available in MACSYMA 

are not implemented on SMP. But users can get the required results by 

implementing those with existing commands. The user's manual for SMP is not 

easy to understand for implementing and using all the capabilities of SMP. 

The people at Inference Corporation are willing to help and guide the users 

to overcome the shortcomings of the manual. The investigators found that 

one of the routines in SMP library, namely FELIP had some bug in it and it 

was reported to Inference Corporation. We believe, it is corrected now. 

Richard Fateman of the University of California at Berkeley, performed 

various experiments with SMP and published various SMP shortcomings of the 

software in ACM's SIGSAM Bulletin of August 1985. These experiments were 

performed on SMP version 1.3 Inference Corporation updated version 1.5 of 

SMP still has some of the old problems. For example, it is unable to factor 

fi fi 

simple algebraic expression X -a . It exhausts memory even if you start 
fresh. This is not an isolated case. Inference claims that these problems 
are scheduled for work prior to next release. 

The other symbolic softwares like REDUCE and MAPLE is also available 
for VAX-11/780 computers. MAPLE and REDUCE, both do not support graphics. 
MAPLE is good for the university environment, in the sense, it does not take 
up lot of memory, and many users can access the same time. One MACSYMA user 
requires 1.5 megabyte compared to 1 OOK for MAPLE. MAPLE is still under 
development and it has limited capabilities compared to MACSYMA. REDUCE is 
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around for a long time and it is widely used in industries and universities 
in U.S. and Europe. Recently, it was announced that REDUCE is now available 
for CDC-CY8ER-170 series. 

Of all different kinds of symbolic manipulation software, according to 
these investigators, MACSYMA has the best capabilities. But the support 
needed for maintenance and development for users of MACSYMA is not available 
from Symbolics Inc. It seems they are more interested in implementing 
MACSYMA on single user machines, rather than supporting MACSYMA software for 
multiuser machines like VAX-11 series. In these respects, SMP software is 
much desirable due to its capabilities and support given by Inference 
Corporation. 

The investigators have performed comparison of execution times of some 
commands in MACSYMA and SMP. 

GENTRAN : 

A software package implementated in RLISP is now available for code 
generation and translation and it runs under REDUCE. It generates complete 
numerical programs directly from REDUCE by transforming REDUCE prefix forms 
into formatted FORTRAN, RATFOR or C code. Assignment statements, control 
structures, type declarations and subprogram headings can be generated from 
algorithmic specifications and symbolically derived formulas. An expression 
segmentation facility breaks large expressions into subexperssions of 
manageable size, and an interface to an existing code optimizer allows 
sequences of assignment statements to be replaced by their optimized 
equivalents. In addition, special file-handling facilities allow code 
generation to be guided by template files and output to be redirected to one 
or more files. GENTRAN provides the flexibility necessary to handle most 
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code generation application under REDUCE. Wm Lei er and Neil Soiffer, of 
Tektronix Inc. has developed a system which is available for workstations to 
do graphics and pretty pointing for REDUCE. REDUCE was chosen over other 
symbolic algebra systems primarily because of its low cost and availability 
of source code. 
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III. REPRESENTATIONS OF SHAPE 


Introduction '• 

In computer graphics and computer-aided design, the problem of shape 
representation and design is a central issue. We are not only concerned 
with pictures but also with providing a geometrical model or database which 
is used for design, analysis and production. Shape representation raises 
many difficulties. There are no well-established criteria for arriving at 
an appropriate mode of description. Approaches to shape and form 
description are usually highly intuitive and adhoc responses to practical 
needs. 

In the early sixties, Charles Coons of M.I.T. initiated some elegant 
and precise methods to represent shapes. His methods were quickly adopted 
by industry; generalized and improved by many researchers. Coons employed 
parametric forms to describe curves and surfaces to obtain axis-independent 
shape descriptions. This also allows the description of closed curves. His 
methods are also piecewise. No attempt is made to represent an entire curve 
or surface by a single analytic function. Curves and surfaces are 
constructed piecewise, then they are sewn together with specified continuity 
condition. P. Bezier proposed a similar idea. Coons' blending functions 
become the basis functions for Hermite interpolation and Bezier's blending 
functions are the basis functions for Bernstein approximation. 

A mathematical representation of shape can oe obtained in two basic 
ways: design and fitting. In design, the problem is to create (or modify an 
existing) shape which satisfies certain constraints. It usually involves 
making interactive changes in shapes and displaying them in real-time. In 
fitting, the problem is to obtain a mathematical representation for an 
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existing shape which is not mathematically defined but exists as a physical 
model from which data points can be measured. 

Traditionally, there are two approaches to these problems: 
approximation or interpolation. Interpolation means the shape matches the 
given data exactly and approximation means one nearly matches the data. 

Most scientific representation of informations requires approximations at 
some level. The approximation might occur at the level of equations that 
model the physical phenomenon, or at the level of the numerical solution of 
these equations. 

Properties of Shape Representation: 

The requirements of the application and of the computer used to 
generate these shapes (curves and surfaces) suggest that their 
representation should have certain properties. Some of the important 
requirements are listed below. 

1. The representation should allow multiple values as in the case of 
closed curves or surfaces. 

0 1 2 

2. Some smoothness is desired. It might be C , C or C continuity. 

3. The representation must not lend itself to too much computation and 
computer memory requirement. 

4. The shape of an object is independent of its orientation with respect 

to arbitrarily chosen coordinate axes: Hence it is desirable that its 

mathematical representation be axis independent. 



5. 


The method to generate these shapes can be local or global. Local 
method is probably preferred because shape can be generated as soon as 
some data points are available. It is also easier to control the 
designed shape by making some modifications around the trouble area. 

6. The mathematical representation should be variation-diminishing, ft 
should smooth small irregularities of the data points instead of 
amplifying them. 

7. It is desirable that the same curve is generated when the data is 
scanned from left-to-right or right-to-left. Hence the representation 
should have syrrenetry property. 

8. The representation should be versatile so that it provides a variety of 
shapes. 

The Representations: 

A survey has been conducted on some of the commonly used shape 
1 -4 

representations . None of these methods requires derivative informations 
directly. Unlike what appeared in many articles on this topic (shape 
representation), the data used in this study are collected physical and 
chemical data and not artificial. SMP is used exclusively to conduct this 
study. Result of this study is described in the following sections. The 
discussion is based on some of the criteria described in the previous 
section. Note that Bezier method is the only approximation, the rest 
provides interpolation. 

Assume there are N data points. They are represented as (x^ ,yi ) , 
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x 2 ’^2 ^ 


• • • » 


P N ( x N »y n ) where x-j < x 2 < . . . < x 


Method A: Oscillatory Method 

Add two dummy data points Pq and P^ + -| » one before the first data point 
and one after the last point. 

2 

For every i = 1, 2, . .., N, define a quadratic polynomial Q.(x) = a x + 
b x + c on [x^, x i+-|] such that Q-j M passes through P^.-j. P^ , P i+ ^. 

On [x., x i+1 ], blend Q^x), Q^x), Q i+1 (x) to form c^x) by using the 
two blending functions as indicated below. 


c.j(x) = 


x ' Vl 


x i x i + l 


x - X, 


Q i (x) + 


x i+l - x i 


Q i+ , (x) 


These polynomical [c^(x) : i = 1 , 2, . . . , N-l] form a piecewise smooth 

interpolant. 


Method B: A Piecewise Cubic Polynomial Method 

Add a dummy point P Q before the first data point. Add P N+1 after the 
last point . 

For every i = 1, 2, 3, ..., N+l, calculate the slope itk of line segment 


y i ‘ y i-l 

m. 

■ i 

x i - x i-l 

An estimated slope s i associated with each point P..(x.j, y^ ) is then 
determined. 

For every i=l,2,...,N, 
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Define a cubic polynomial function c^(x) on [x^, x. + ^] such that c^(x) 
interpolates P^_^, P.. , P..^ and c.j (x^) = s^ . 

Piece these polynomials [c^(x) : i = 1, 2, .... N-l] to form a smooth 
interpolant over [x^, x^]. 

Method B ' : A Piecewise Cubic .Polynomial Method with Akima's Slope Estimate 

This method is identical to Method B except that the estimated slope s^ 
at point P i (x i , y.. ) is defined as: 

S 1 = m 2 

s 2 = 0.5 (m^ + m 3 ) 
and for every i = 3, N-2 

' m i+2 ' "l+l 1 ”i * 1 m i-l ' m i 1 Vl 

s i = 

|m H2 • "Vl I + 1 m i-l * m 1 1 
which is a weighted average of m.. , nr^. 

S^_] = 0.5 (m^_i + m^) 

S N = m N 

Method C: A Piecewise Cubic Hermite Polynomial 

The slope S i at data point P^x^, y . ) is estimated as in Method B. On 
each interval [x^ , x ^ + ^ ] , a cubic polynomial c ^ ( t ) in parameter t is 
constructed by using blending functions as indicated. 

c i (t) = ( 2t 3 - 3t 2 + 1) P. + ( -2t 3 + 3t 2 ) P. +1 + 

(t 3 - 2t 2 + 1) p: + (t 3 - 1 2 ) p^ +r 
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where 0 £ t £ 1 . 

These c^'s , i = 1, 2, . .., N-l give a smooth interpolant. 

Method D: McLaughlin's Method 

Each data point with the exception of the last point is classified as 
"regular"; or "irregular". If a point P^x^, y^ ) is irregular, a line 
segment is drawn from this point to the next point P i+1 (x i+ ^, y i+1 )- If a 
point is regular, a break point B(x', y 1 ) will be inserted between P. and 
P^ + .j, then a parabola passing through P.. and P. + ^ will be drawn. 

Method E: Piecewise Bezier Method 

Four control points will be used to construct a Bezier curve. 

3 

P(t) = £ (jj)t k (l-t) 3_k P. +|< 
k=o 

where 0 <_ t <_ 1 . 

Piece these curves together to get an approximation. 

Data: 

The data used to generate Figures 1 to 20 are given in Tables 1-5. 
Water samples are collected from the same station at a depth of one meter. 
Some physical measurements are taken. The samples are also analyzed for 
some of their chemical contents. 

Variable x represents time in the unit of weeks, y represents the 
measurement of a physical or chemical variable at low tide or high tide. 
Table 1 shows Dissolved Silica content (in mg/1 ) at low tide (yl ) and high 
tide (y2). Similarly, Tables 2-5 show the measurements of Salinity 
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Conductivity, pH, Total Kjeldahl Nitrogen and Total Phosphate in appropriate 
units at low tide and high tide. 


X 


yi 

y2 

1 


1.4 

1.3 

2 


1.8 

0.3 

5 


1.7 

1.6 

6 


1.6 

1.7 

7 


1.6 

1.5 

9 


0.3 

0.1 

11.5 


0.1 

0.8 

14 


0.2 

1.8 

17 


0.9 

2.2 

22 


0.4 

0.8 

Table 1. 

Dissolved Silica 

in mg/1 


X 

: time in weeks 


yi 

: low tide 

data 


y2 

: high tide 

data 

X 


yi 

y2 

1 


20.5 

21.1 

2 


23.1 

23.2 

5 


25.2 

26.0 

6 


27.6 

25.5 

7 


19.8 

19.7 

9 


20.9 

29.1 

11.5 


26.4 

23.1 

14 


23.8 

21.4 

17 


25.0 

21.1 

22 


23.9 

24.6 

Table 2. 

Salinity Conductivity in 


x : time in weeks 

yl : low tide data 

y2 : high tide data 
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X 


yi 


y 2 

1 


8.45 


8.28 

2 


7.98 


8.45 

5 


8.62 


8.31 

6 


7.87 


8.16 

7 


8.06 


7.81 

9 


8.46 


8.23 

11.5 


7.77 


8.21 

14 


7.68 


8.10 

17 


8.12 


8.23 

22 


8.25 


7.98 

Table 3. 

PH 

in standard 

unit 


X 

: time 

in 

weeks 


yi 

: low tide data 


y2 

: high 

tide data 

X 


yi 


y2 

1 


0.351 


0.909 

2 


0.158 


0.342 

5 


0.698 


1.173 

6 


0.868 


0.884 

7 


0.765 


0.571 

9 


0.500 


0.621 

11.5 


1.320 


1.084 

14 


0.552 


0.705 

17 


1.912 


1.111 

22 


1 .088 


0.533 


Table 4. Total Kjeldahl Nitrogen in mg/1 
x : time in weeks 

yl : low tide data 

y2 : high tide data 


X 

yi 

y2 

1 

0.134 

0.182 

2 

0.059 

0.157 

5 

0.065 

0.042 

6 

0.052 

0.053 

7 

0.080 

0.100 

9 

0.042 

0.060 

11.5 

0.190 

0.192 

14 

0.018 

0.034 

17 

0.157 

0.171 

22 

0.041 

0.000 

Table 5. 

Total Phosphate in mg/1 


x : 

time in weeks 


yi : 

low tide data 


y2 : 

high tide data 
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Conclusion : 


The first three methods (A, B, B‘) desbribed in Section III define 
polynomials on [x^ , x. -|] using explicit form 
y = f(x). 

Hence none of them permits multiple-valued shapes. Cubic Hermite form, 
Bezier curve and McLaughlin's method provide curves of single-valued 
functions as well as multiple-valued curves. 

The oscul atory method (A) matches first order derivatives at the two 
end points of interval [x i , x i+1 J as well as the data points P^x^ y. ) and 
P-j + -| ( x i + -| > y-j+i^* ^ orc * er This can be seen from the following 

direct calculation. From Section III, 


c i (x) 


x x i + l 


x i " x i+l 


x - x. 


P-j (x) + 


<W x) on [x,, x 1+ ,] 


x i+l ' x i 


(x) 


q-i (x) 


x - X 


i+1 


X i " x i+l x i ' x i+l 


q i+ ,(x) 
q.j(x) + + 


x - 


x i+r x i x i+l " x i 


q- tl (x) 


c i (x i> 




x i " x i+l 


+ q:(x.) + 




x i+1 ‘ x i 


q-i (x i) - WV 


x i " x i+l 


y, - y ( 


+ «i (x i ) 


x i * x i+i + 


W • 
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q i+l (x i+l } 


Similarly, 


c i (x^+1 ) 


Methods B, B' define curves over [x^, x^ + -|] in such a way so that only , 
P ^ + 1 and the first order derivative at x^ are matched. As a result, they 
only provide curves. Method C generates piecewise cubic Hermite 
polynomials which are continuous. McLaughlin's method is continuous 
since derivatives are not considered directly by this method at all. Bezier 
method can supply continuity of any order provided the data points satisfy 
certain conditions. In this investigation, the data poins are collected 
data and will not satisfy those conditions. Hence the Bezier curves are 
only C° continuous. 

All methods described in Section III are local methods. Bezier method 
is actually a global method. But in this investigation, a piecewise Bezier 
curve is used instead of a single Bezier curve therefore it is considered as 
a local method. 

Bezier curves are symmetric due to the symmetry of the Bezier blending 

functions. Take four consecutive data points P., P., n , P.,„, P-.-,, they 

r l l+l i+2 i+3 J 

generate 

3 

P(t) = £ (j|) t k d-t )3Hc p. +lc . 0<t<l 

k=o 

The points P i+3 , P i+2 , P i+i» P i in this order (from right-to-left) generate 
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R(t) = D (^ } tk n-t) 3_k p i+3 _ k . o< t < l 

k=o 


Take any to such that 0 <_ t 0 <_ 1 


3 

p(t o> - s ( k> C "V k P i + k 

k=o 

3 

Rd-t 0 ) . 2 (*) n-t 0 ) k t o 3 ‘ k P i+3 . k . 

k=o 


Let j = 3-k, then 


0 

»n-t 0 ) - £ ( 3 -j> <i-t 0 ) 3 ' J V p 1tJ 

j=3 

0 

L ( 3 ) t 3 (1-t ) 3_J P. . 

V 0 v o' 1+J 

j=3 


= P(t 0 ). 

McLaughlin's method is not symmetric. The determination of the break point 
B depends on some "previous" points. While moving from right-to-left along 
the x-axis, these "previous" points are not the same as those associated 
with moving from left-to-right. In general, different parabolas will be 
constructed between P. and P^. However, the shape of the generated' curve 
is consistent in the sense that if is determined as regular and a 
parabola is found between P.. and P i+1 , when traversing from right-to-left, 
P-j +1 will be determined as regular and a parabola (a different one) will be 
constructed between P. + i and P. . 


Cubic Hermite form is also symmetric. From P^ to P^ + p one has curve 


19 


represented as: 

P(t) = (2t 3 - 3t 2 + 1) P(0) + (- 2t 3 + 3t 2 ) P(l) + 

(t 3 - 2t 2 + t) P ' (0) + (t 3 - t 2 ) P' (1) 0 < t < 1. 

From P i+1 to P^, one will have 

R(t) = (2t 3 - 3t 2 + 1) P(l) + (- 2t 3 + 3t 2 ) P(0) + 

(t 3 - 2t 2 + t) (- P ' ( 1 ) ) + (t 3 - t 2 ) (- P ’ (0) ) 0 < t < 1. 

For every t in [0,1] 

R(l-t) = [2(1 -t ) 3 - 3(1 -t ) 2 + 1] P(l) + [- 2(1 -t ) 3 + 3(l-t) 2 ]* 

P(0) + [(1-t) 3 - 2 ( 1 -t ) 2 + (1-t)] (- P 1 ( 1 ) ) + 

[(1-t) 3 - (1-t) 2 ] (- P 1 (0) ) 

= (2 - 6t + 6t 2 - 2t 3 - 3 + 6t - 3t 2 + 1) P ( 1 ) + 

(- 2 + 6t - 6t 2 + 2t 3 + 3 - 6t + 3t 2 ) P(0) + 

(1 - 3t + 3t 2 - t 3 - 2 + 4t - 2t 2 + 1 - t) (-P'(l)) + 

(1 - 3t + 3t 2 - t 3 - 1 + 2t - t 2 ) (-P'(O)) 

= (- 2t 3 + 3t 2 ) P(l) + (2t 3 - 3t 2 + 1) P(0) + 

(t 3 - t 2 ) P'(l) + (t 3 - 2t 2 + t) P ‘ (0) 

= P(t). 

The remaining methods do not have the symmetry property. 

A Bezier curve is independent of the coordinate system used to measure 
the locations of control points. The same is true for the McLaughlin's 
method because both use parametric form and the methods involve only control 
points (and also beak point for the latter). Cubic Hermite form is not axis 
independent despite parametric form is also used. It depends on first order 
derivatives as well as data points. These first order derivatives do not 
undergo the same transformation as the axes. The other methods are also 
axis dependent. 

The piecewise Bezier method is only C^. Curves generated by this 
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method in general contain corners and/or cusps. Nonetheless the method 
provides very interesting curves. As can be expected, McLaughlin's method 
generates too many line segments along the process. Some curves are- made up 
of line segments only. However, when some line segments are mixed with a 
few parabolas, the results are quite interesting and pleasing to the eyes. 
These two methods could be valuable to those applications that call for 
corners or cusps. The cubic Hermite method does not produce satisfactory 
results. The curves have too many twists and turns. Method A gives the 
smoothest curves. On a few occasions, the curve does overshoot data points. 
Method B, B' give very similar results. Curves are relatively smooth, 
though not as smooth as those of A. They also produce some corners. The 
overshoot parts are not as noticeable as those of A. Some curves produced 
by B' appear to be a bit tighter in some areas. 

The approximate amount of computer time in the unit of seconds needed 
to generate the curves of Figures 1-20 are given in the following table. 
These curves are generated by SMP version 1.5.0 running on a VAX 11/780 and 
a DEC GIGI is used as graphic display. 
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METHOD 



A 

B 

B 1 

C 

D 

E 

DS 

12.03 

11.15 

11.95 

33.12 

22.35 

22.88 

SC 

12.05 

10.98 

11.78 

33.12 

15.83 

22.77 

pH 

12.10 

11.08 

11.71 

33.13 

18.35 

22.85 

TKN 

12.00 

10.97 

11.65 

33.20 

19.85 

22.73 

TP04 

11.83 

10.73 

11.51 

33.25 

13.27 

22.75 

DS 

11.97 

11.07 

11.68 

33.33 

25.20 

22.68 

SC 

12.05 

11.02 

11.68 

33.15 

26.57 

22.90 

pH 

11.90 

11.28 

11.57 

33.22 

22.78 

22.90 

TKN 

11.91 

10.98 

11.62 

33.07 

18.43 

22.65 

TP04 

11.97 

10.80 

11.70 

33.33 

15.85 

22.02 


Table 6. Approximate time needed in sec to generate curve 
Abbreviations: DS - Dissolved Silica 

SC - Salinity Conductivity 

TKN - Total Kjeldahl Nitrogen 
TP04 - Total Phosphate 
1 - low tide 

h - high tide 
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Figure 2. Dissolved Silica at low tide 
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Figure 18. Total Kjeldahl Nitrogen at high tide 
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SMP Code: 


/* This program uses the osculatory method on page 11 of 
Brodlie's Mathematical Methods in Computer Graphics 
and Design to draw a smooth curve that passes through 
10 given data points. 

An additional point is added to each end of the data. 
It is the mirror image of the first point (or the 
last point) of the data. 


A quadratic polynomial q( i , x ) =ax A 2 + bx + c defined on 
[ x ( i ) ,x(i+1 )] is determined. This polynomial passes 
through points P(i-l), P(i), P(i+l) and is given by 
_y( i-1 ) 1 [x( i-1 ) A 2 x( i-1 ) ll [a. 

.y(i) . = ,x(i) A 2 x ( i ) 1 . .b. 

_y(i+l)J [x(i+l) A 2 x(i+l) 1J [c. 

On [ x ( i ) , x ( i + 1 ) ] , blend q(i,x), q(i+1,x) to get c(i,x) 

by using 

x - x( i + 1 ) x-x(i) 

c(i,x) = q ( i » x ) + q(i + 1,x) 

x(i) - x(i+l) x(i+l) - x(i) 

Join these {c(i,x) : i = 2,3, • • ,10} into a piecewise 
smooth curve */ 

x : {-1,1,2,5,6,7,9,11.5,14,17,22,27! 

D°[i,2, 1 , \ 

tl : Uyli- 1 1 } , jy[i] ! , |y[i+1 1 } h \ 

1 2 : x[i-1 ] A 2 , x[ i-1 ] , 1 ] , { x[i] A 2 , x[ i ] , 1 | , { x[ i + 1 ] A 2 , x[i + 1 ] , 1 } | ; \ 


tl : { yLi-i J} » tyLiJI tiyLi+i Jti; \ r , , , r 

1 2 : | { x[i-1 ] A 2 , x[ i-1 ] , 1 } , { x[i] *2 , x[ i ] , 1 | , { x[ i + 1 ] A 2 , x[ 

coef : Mdiv[t1,t2j: \ 

q[i,$u] : c o e f [ 1 , 1 ] ($u) A 2 + c o e f [ 2 ,■ 1 ] $u + coef[3,l]] 


Do[j,2,10 \ 

c[j,$u] : ( ( $u - x 

( ( $u - x 


. j + 1 ])/(x[jl - x[ j + 1 ])) q[ j,$ul + 
.j])/(x[j+1J - x[j])) q[j+1,$uj] 


f t $u] : c [ 2 , $u] 
f _$u_=$u>2 i : c 3,$u_ 

f,$u_=$u>5. : c_4,$u_ 

f.$u_=Su>6_ : c_5,$u_ 

f > $u_=$u>7. : c > 6,$u_ 

f,$u_=$u>9. : c_7,$u_ 

f . $u_=$u> 1 . 5 ] : c[8,$u] 

f i $u_=$u>14. : c_9,$u] 

f F $u_=$u> 1 7 1 : cF 1 0 , $u] 

_Graph[Smp]fl ] : 0 
_Graph1 Smpj[l]:0 
Graph[f[t],.t,1 ,22, , , , { { 0 , 22 } } ] 
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I* Method B */ 

x : {-1 f 1 ,2,5,6,7,9,1 1 .5, -14, 17,22,27} 


/* Estimate slope at x(i) */ 

Do [ i , 2 , 1 2 , \ 

m[i] : ( y [ i ] - y [ i- 1 ] ) / ( x [ i ] 
Do[j,2,11, \ r n 

s[j] : 0.5 (m[ j ] + m[ j+1 ] ) ] 


x[i-l])] 


/* Determine cubic polynomials c(i,x) */ 
Do [ k , 2 , 1 0 , \ 

tl : Jy[k-1 1 } , |y[k] } , {y[k+1 1 } , { s[kl }} 

1 2 : (x _k-1 ] *5 , xrk-l] A 2, x.k-l],l] f 

{ x k ] A 3 ,xk] A 2 , x k ] ,l' 

| x[ k + 1 ] A 3 , x[k + ll A 2, x[k + l], 1 
i3*x[kj A 2, 2*x[k] ,1 ,0 

coef : Mdiv[t1,t2] ; \ 
c[k,$u] : coeff 1 , 1 _ ($u) A 3 + coef[2,l] 

coef[3,1. $u + coef[4,1j] 


\ 


/* 

f 

f 

f 

f 

f 

f 

f 

f 

f 


piece c ( i , x ) 
$uj : c[ 
$u_=$u> 2 
$u_=$u> 5 
$ u_= $ u > 6 
$u_=$u>7 
$ u_= $ u > 9 
$ u_= $ u > 1 . . . 

$u =$u>1 4] 


together */ 


C , MJU J 
: C 

3 , $u 

: C 

4, $u 

: c 

. 5 , $u_ 

: 0 

6 , $u 

: c 

>,$ u] 


^ l 9 ^ 

: [ 9 , $u] 


_$u_=$u>17J : c[lO,$u] 

_Graph [ Snip ] [ 1 ] : 0 

Graph [ f[ t J , t , 1 ,22, , , , { { 0 , 22 } }] 


; \ 

\ 

\ 

; \ 

( $u ) A 2 + \ 
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/* Method B' */ 


x : !-1 , 1 , 2, 5,6,7, 9, 1 1 .5, 14, 17, 22,27} 

/* Calculate Akima's slope estimate at x(i) */ 

Do [ i , 2 , 1 2 , \ 

mfil : (y[i] - y [ i- 1 ] ) / ( x [ i ] - x[ i - 1 ] ) ] 
sf2l : m[3] 

s[3J : 0.5 (m[ 3 ] + m[4]) 

Do[j,4.9. \ 

s [ J ] : (Abs[m[ j + 2 ] - m[j+l]] m[j] + Abs[m[j-l] - 


1 [ 1 0 ] : 

s L 1 1 : 


( Abs[m[ j + 2 ] - m[j+l]] m.j] + Abs[m[j-l] - m[j]l * \ 
®[j+l])/ (Abs[m[j+2j - m[j+l]] + Abs[m[j-l] - m[j]])] 
: 0.5 (m[ 1 0 ] + m[l 1 j) 


/* Define cubic polynomials c(i,x) */ 

Do[ k , 2 , 1 0 , \ 

tl : |y[k-1 } , }y[k] } , jy[k + 1 1 } , { s[k] } } ; \ 

1 2 : x[k-l] A 3, x.k-1J A 2, x.k-l],l], \ 

x[k] *3 , x ,k] A 2 , x}k] , 1 , \ 

x[k+lJ A 3, x[k + l] A 2, x[k + l], 1 . \ 

3*x[ k J A 2 , 2*x[ k ,1 , 0 j ; \ 

coef : Mdiv[t1,t2j ; \ 

c[k,$u] : coef f 1 , 1 1 ($u) A 3 + coef[2,l] ($u) A 2 + \ 

coef[3,1J $u + coef[4,l]] 

/* Piece c(i,x) together */ 
f [$u] : c[2 , $u] 

f f $ u_= $ u > 2 _ : c r 3,$u_ 

f[$u_=$ u >5. : c_4,$u_ 

f $u_=$u>6. : c^.Su, 

f $u_=$u>7. : c_6,$u_ 

f $u_=$u>9. : c_7,$u. 

f[$u_=$u>1 1 . 5 ] : c [ 8 , $u] 
f [ $ u_= $ u > 1 4 1 : c f 9 , $ u ] 
f [ $u_=$u> 1 7 1 : c [ 1 0 , $u ] 

_Graph[ Smp] [ 1 ] : 0 
Graph[ f [ t J , t , 1 , 22 ] 
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/* Cubic hermite form */ 


x : {-1,1,2,5,6,7,9,11.5,14,17,22,27} 


/* Estimate first order derivative at x(i) */ 
Do [ i , 2 1 2 , \ 

mfi] : ( y [ i ] - y [ i- 1 ] ) / ( x [ i ] - x[i-l])] 

Do[ j , 2 , 1 1 , \ 

s[j] : 0.5 ( m[ j ] + m[ j + 1 ] ) ] 


P * { I 

Do [ k , 1 , 1 1 , 


q : 1 I 

Do [ k , 1 , 1 1 , 


P 

q 


Cat[p, { {x[k] ,y[k]}} ]] 
Cat[q, { { x[k] , s[k] } } ] ] 


/* define the cubic polynomials */ 

Do [ k , 2 , 1 0 , \ 

f[k,$u] : (2 ($u) A 5 - 3 ( $u ) A 2 + 1) p[k] + \ 
(-2 ( $u ) *3 + 3 ($u) A 2) p[k+ll + \ 
(($u) A 3 - 2 ( $u ) A 2 + $u) q[k J + \ 

( ($u) A 3 - ( $u ) A 2 ) q[ k + 1 ]] 

Graph[{f[2,tl ,f[3,tl ,f[4,t],f[5,t],f[6,t],f[7,t] , \ 
f[8,t],f[9,t J,f[lO,t] } ,t,0,1 ] 
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/* McLaughlin's Method */ 
delta[$i] : p[$i] - p [ $ i — 1 ] 


/* Determine the regularity of the first point */ 


regularl[$p] : [ a 

b 


$p[l 

$PL3J 


a . b > 0 


_ « ' 
- $p[2. 


\ 

\ 


/* Determine the regularity of the last point */ 


regularN[$p , $n] : [ 


a 
b 

a . b > 0 


$pJ_$n-2 j - $ pj_ $n- 1 


$pL $n J - $p|_ $n- 1 J 


\ 

\ 


/* Determine the regularity of the in between points */ 


regular [ $p , $i 

[ d 


\ 


d 1 


Detf { {Elem[delta[$i 
]Elem[delta[$i] , 
Det[ { {Elem[delta[$i 


,{l)], Elem[delta[$i+2] , [l }]} 
2}].Elem[delta[$i + 2], { 2 }]} 


J . 

,d 1, 


El em 
El em 


delta[$i+1 
delta[ $i 


{Elem[delta[$i], { 2 } ] , Elem[delta[$i + 1 1 



Elem[delta[$il, 

d2 : De t [ { { Elem[ delta[ $i + 2 ] , { 1 } ] +Elem[ delt a[ $i + 1 j , 

Elem[delta[$i + 2],[l } J 
{Elem[delta[$i+2], {2}]+Elem[delta_$i+1 1, {2 

Elem[delta $i + 2J, {2 

del[$i] : d ; \ 

dels[$i] : dl ; \ 

And[d ~= 0, Sign[d] = Sign[dl], Sign[d] = Sign[d2]]] 


norm[$a,$b] : ($a[l] - $b[l]) A 2 + ($a[2] - $b[2]) A 2 

/* the set of data points P has been batched in already */ 


q : 0.5 
n : 1 0 



rl : Li s t [ regularl [ p ] ] 
r : Cat [ r , rl ] 
rl : Li s t [ regular [ p , 2 ] ] 
r : Cat [ r , rl ] 


t[l] : If[ 



(p[3] ; P[2]> 


- P 

" P 


) = \ 


).(P 1 
j).(p[2_ 

e2 : norm[a,p[2]] = norm ta, P [l]] 
Sol [ { el , e2 } , sj ; \ 


ff] 


) ; \ 


\ 
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i t [ 1 ] : P [ 2 ] + s*( p[?] - P [ 2 ] ) ; \ 


a 

b 

fl 


p[i 

PL 2 


+ q*(it[l J - p[ 1 J ) ; \ 


( b - p 
( b - p 


h*( 


a - 


P [ 2 ] ) ; \ 


) • (p[l . - P 

) . ( p [ 2 J - p 


'2 ' 

'l ' 


) = \ 
) ; \ 


r L J / • \ .r L J FL ' J / 9 

f2 : norm[b , p[2 J J = norm[ b , p[ i ] ] 
So 1 [ { f 1 , f 2 } , h J 


[ 


p[2] 
a : 
e 1 : 

e2 : 
So 1 [ 
it[l 


h*( 


p[2] 


a - 


J ; \ 

P 21)] 
( P L 3 J — 


P [ 2 ] ) 


\ 

) - \ 
) ;.\ 


p[i ]] 


i t [ 1 


7 - r n - L ^ J r _ L J / 9 

(a - p..2. ) .(p[l 1 - p f 2 1 
(a - p[l 1) . ( p [ 2 J - p[l J 
normL a , p[2 J J = norm[a,p[l]J 
el , e2} , s J ; \ 

• p[ 2 J + s*(p[5] - p [ 2 ] ) 
]]] , \ 


Do[ j » 2 , n-2 , \ 

rl : If[ j<n-2, regularf p , j + 1 _ 

regularN[ p , n] _ 

r2 : Li s t [ rl ] ; 
r : Cat [ r , r2 ] ; \ 

If[ r[j]=1,. . \ 


[ If 


r[j + l]= 1 , \ 

s : 1 + dels[ j] / del 

i t [ j ] : p[ j+ 2 J + s*(p 


\ 

\ 

\ 


^ , 
j + 1 


[j] + q* ( i t [ j ] - p[j1) 


- p[ j + 2 ] ) 


PL J + 1 ] + s*(a - p[ j + 1 ] ) 


[ 


a 
b 

c : p[j] + h*(t[ j-1 ] - p[ j]) 
Sol[b=c , h] ; 
t[jJ : p[j] + h*(t[j-l] - p[j])], 
1 + dels [ j ] / del[j] 


t[n-1 ] 


t[ j ] 

If C 



p[j+2J + s*(p[j+l] - p[j + 2 ])]]], \ 


n_1 ] = 1 i \ n r 

a : p[n-l] + s*(t[n-2 


el 


(a - 
(a - 


EV 


).(p 
) • (p 


n 


- P 


n-1 1) 


J - PL n “1 . ) = 
n-1 ] - p[n] ) 


e2 : norm[ a , p[n-1 ] ] = norm[a, p[n] ] 
So 1 [ | e 1 , e 2 } , s 


E:::]]* 


PL n ~ 

p 


S*( t[n-2] - p[ n-1 ] ) ] 


\ 

\ 

\ 

\ 

\ 

\ 


\ 


Do[ j , 1 , n-1 

f[j,$u] : \ 

If[r[j] = 1, $u A 2 p[j + l] + 2 $u ( 1 — $ u ) t[j] + \ 
( 1 -$u ) A 2 p 


>L j 


p[ j] + (pJi + i ] - p [ j ] ) ] ] 

Graph [ ff fl ,$ul ,f[2,$u] , f[3,$ul ,f[4, $u] , f [ 5 , $u ] , f [ 6 , $u ] , \ 
f[7,$u] ,f[8,$uj,f[9,$u]} ,$u,0,1 ] 


\ 
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/* Bezier curve using 4 control points */ 
x : {1 ,2, 5, 6, 7, 9,1 1 -5,14,17,22} 


x[i],y[i]}}]] 


Do[i, 1 ,7,3, \ 

f [ i , $u] : Sum[Corab[3,k] ($u) A k (l - $u)*(3-k) p[i+k], 


q : Map[ Pt , p 

Graph [ ] f[l , t J , f [ 4 , t ] , f [ 7 , t ] } , t , 0 , 1 ] 


{k, 0, 3 } ] ] 
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IV. KINEMATICS OF ROBOT ARM 


Introduction : 

The purpose of this research work is to study the kinematics of a 
single link robot arm. In the solution of equations of motion of this 
system, the symbolic manipulation programs MACSYMA and SMP have been used. 
Also, some IMSL subroutines are also utilized. 

The Lagrangian of the system is set-up in a fixed frame of reference. 
Hamiltonian is calculated. From this, Hami Iton-Jacobi equation is obtained. 
Solution of this equation gives the time evolution of the system. 

Definition of symbols: 

Before setting-up the necessary equations, the follwoing list of 
symbols is defined. Note that subscript 1 denotes the first rod and 
subscript 2 denotes the second rod. Subscripts to I denote the appropriate 
axes about which the moment of inertia of a particular rod is defined. See 
Figures 1 and 2. 

L = Lagrangian 
T = Kinetic Energy 
V = potential energy 
I = moment of inertia 
<t> = angle of rotation 
m = mass 

g = gravitational acceleration 
h = height 
w = angular speed 

1 = distance of CM from axis of rotation 
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Lagrange's Equations: 

The Lagrangian of the first rod is 

h = T 1 - V 1 


1 T 

2 I 1Z ] Z 1 ♦l 


m i g h ] 


Kinetic energy of second rod in the X^Y^Z^ system is 


T = — T •? 

'2 2 1 2Y 2 Y 2 ^ 


Now the vector is 


1 — d? y 

2 ? 2 y 2 


In X^Y-jZ^ system, the vector | 2 is given by 


(d<t> 2 /dt) x ^ Y ^ z ^ = (d<t, 2 /dt) X 2 Y 2 Z 2 + W 1 x *2 


This can be shown to give 


(<J 2 ) 

2 X^Y ]Z] 


*2 •? ? 
*2 + *1 *2 


Therefore, 

T 2 in Y-j Zi system is 


(TJ 


1 1 i 2 + 1 1 


•2 , 2 


2 / X ] Y 1 Z 1 ~ 2 a 2Y 2 Y 2 . T 2 "2 A 2Y 0 Y„ *2 


2 2 
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and 


(V 2 ) X ] Y 1 Z 1 " m 2 9h 2 ' 


Kinetic energy of motion of CM of rod 2 about Z-j axis is 


(T 2 ) Z ] ' 2 m 2 1 2 *1 * 


From above, the Lagrangian of second rod in the XjY-^ system is 


1 2 *2 1 *2 

2 m 2 1 2 *1 + 2 I 2Y 2 Y 2 *2 


+ 1 T :2 2 

2 X 2Y 2 Y 2 *1 *2 


m,gh 


2 a 2 


Total Lagrangian can be shown to be 

L - L i + 4 


= I 


2 .2 


2 i *] + 2 ^2 ^2 + 2 ^2 ^1 ^ 


where 


*1 ' 2 (I 1Z ] Z 1 + m 2 1 2 } 


I = — I 
2 2 2Y 2 Y 2 


V = g(m^ h-j + m 2 h 2 ) 
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Lagrange's equations of motion are 


<f> 

1 


<t> 

2 


2 *2 

(a + 4>2 ) 


4> 

2 


2 

1 


wi th 


Hamiltonian i-ormulation: 

Let the generalized coordinates be q = ( ^ , <t> 2 ). Generalized momenta are 
vector p. 


Now, 

P-, = 3L/ aq-j = ^ ^ + i 2 ^ ^ 

P 2 = 3L/aq* 2 = l 2 <J 2 . 

Hamiltonian H is 
2 

H = £ q i P i - L 
i =1 11 

1 ■, *2 1 *2 1 *2 2 

= 2 ^1 ^1 + 2 ^2 ^2 + 2 ^2 ^1 *2 + V 

= p^/(2( I + I 2 <d|) ) + p|/(2I 2 ) + V. 
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Hamilton's equations of motion are 

A 

<t>1 — 9 H / 3 p i — p i / ( 1 1 Ip $2 ^ 

$ 2 = aH/ ap 2 = p 2 /I 2 

p-| = - aH/ aq -j = 0 

P 2 = - 3H/ 9q 2 = I 2 i> 2 P^/ ( 1 -j + I 2 4> 2 ) 2 . 

Note that p 1 = constant. 

If one assumes the system to be non-dissipative, H = E = constant. In 
this case. 


p = + a*(1 - — - — p-) ^ , 
d ~ i + q 

remembering that p-j = constant and defining 
C = B/A > 0. 

with 

A = 2 (E -V ) 


and 

B = y 2 (= p 2 ) > 0. 

This leads to the interesting phenomena of bifurcation. When 0 <_ C < 1, 
then - » £<t^ _< “ . This corresponds to rotation. When C > 1, ^ i ' (C-l)^. 
This corresponds to libration. Thus C = 1 is a bifurcation point. 
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Hamilton-Jacobi Equation: 

Redefine the momenta and coordinates as 


PT = P-j 


p 2 = p 2 /l 2 c 




From now on, p^, p 2 and $ 2 will denote p^, p 2 and 4> 2 - With this 
redefinition. 




+ ‘ + V. 


2(1 + <t> 2 ) 2 

Again, assuming the system to be conservative, Hamilton-Jacobi equation will 
be given by 


p 2 

+ — + V 


2(1+ <t> 2 ) 2 


with 


3W 
3 4 > ^ 


From Hamilton-Jacobi theory, one finds that 


W = Y<f \ + f (A - — y-)* d<t>p . 

1 J (1 + * 2 ) 2 
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This leads to 


Q 1 ~ t + 6 1 

= 3W/3E 



6 1 is to be determined from initial condition on <t> 2 . 


Simi 1 arly. 


Q 2 - 6 2 


= 3 W/ 3 Y 


which yields 


'1 = B 2 + * 


/ 


d<t>. 


(1+ *\) [A - — 

1+ 4>o 


Again, 6^ has to be determined from initial condition on 4>^ 


In terms of C = B/A, one obtains 


t + B, = -V 


and 


■if 


d<t> , 


[1 - : " C 2 3^ 

( 1 + 2) 


= + 


Y 


f— 

J (1 + 


d<t> 


2 


2 (1 + $ 2 ) [ 1 - — < £— ]£ 

1 + 4>2 
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Evaluation of integrals: 

If one makes a change of variable from 4 > ^ to e with <t >2 = cot e, the 
first integral reduces to 

-F ( e , K ) + E (e,K) + a cot e 

where 



a = [1 - C sin^e]^ 
and 

F = Elliptical integral of first kind 
E = Elliptical integral of second kind. 

The second integral can similarly be evaluated to -F(e,k). 

This completes the solution of the equation of motion of system. With 
the evaluation of these integrals, one can obtain tU^). This can be 
inverted to give ^(t). Once ^(t) is known, <t>-|(t) can be obtained from 
second integral. 

Use of MACSYMA: 

The MACSYMA program which obtains the Hami 1 ton- Jacobi equation is shown 
on the following pages. In Figure 3, is shown as a function of time t. 

To evaluate the integrals, first SMP was used. But in doing so, some 
serious errors in the elliptical functions evaluation programs of SMP were 
discovered. Inference Corporation was informed this. SMP has now corrected 
these errors. Since SMP could not be used to evaluate the elliptical 
functions, the subroutines MMLINF, MMLIND and MMLINJ of IMSL were utilized. 


57 



17 11 


ORIGINAL PAGE IS 
PF POOR quality; 

1 J ;; e p ; .. J 1 


fc 


■ . 3 i 1 J :: o 


: -• s n d cd 12, £ a 1 > a23 , Eh3 i Cal » u 2 :> 1 - 

l 2 ) / < . ' r ; i 2 i r i , - 7 


n i ' f ( 


- » c. * .■ ? 


( d i } 


• ■ ca t'll , t > ; 


.-■2 2 Cal » o 2 r r 1 - r-2 2 C t 2 ) J 

1 1 £ ■ t Cl 2 < t > r r- 1 ( t > ♦ p 2 < t ) 'J 


«. « f •- 


GCt 1 


i r : ) 


d c. : 

ct 


<dd) 


v - * • \ u 4. i, j r 


{ ^ r.* * •/ - 

• - * ' •- ^ ; * 


ac? : 


dt 


( d 2 ) 


- / — *. A. 


d t 


• ■- •' - " 2 + 1 / 2 :r i 2 d c 2 ' 2 1 / 2 V i 2 d d c 1 - 2 l g 2 - 2 - v ? 


(22 > 


d a 2 2 t d a 1 2 2 

2 2 ( > 12 ( ) u 2 ii 

■~u dt 

+ + -- 


( ) 

' dt ‘ 


<c-?> di a (ciff ( I r d •?, l ) , t )-diff( 1 , &. l ) — 0 ? 


dt 


f - i A 


- d 1 v o ( *: , a : f f ( d a irt) ) f 


58 


( diO ) 


(e 11) 


ORIGINAL PAGE IS , 

OF POOR QUALITY. •'■ 5 ' 'L' 

- , .. c d : « i . J .?;X' : Z 0 & 


12 ' a 2 + u 


d i f f ( d i f f ( 1 f do2 ) » t ) — cf i ff ( } , q2 ) =0 > 


(dll) 


d q 2 . 


dal 2 


i 2 - 12 ■•< — ) ’ *: 


* dt •; 


(cl2) • soIve(Zrtiiff <da2vt) : ) ; 


( d 1 2 ) 


d Ct 2 • rj C( l 2 ; Tr" ' 

C *.<••--«) *22 

. 2 dt 

"dt -'••-••••• — 


, . . . 1*V 

. . r *. ** • - ! , |<" ... . ' 


li.. , 


( cl3) F'l"il*dal-K2*dal*a2"2? 


C d 1 3 ) 


(cl-)) so 1 ve ('/. * dal ) ; 


'( d 1 -i ) 


(clS) F-2=i2#dc2 i 


~(dl5 ) 


(eld) solve (X r da2 ) ! 


“(did) * 


. { cl7) >srt(dl3»2) » 


dal -2 •• dol 

:>1 * i2 — - a2 - + :.ii - 

'■ "dt- ~ •; — dt 


dal pi ; • 

C W - 

' dt ~ 2" : r 

i2 a 2 fTf ;- 


da2 

T '2 -" 12 ' '--~y 

■ dt .... • i 


da2 ' 

•• *"C~~ ~ “'"'' 2 ' "’r ; 1 
dt • i2 ■■ r 


■ . • - ..v.);. ' : v’ id 

v v. • >• ; 




. • '• • . ... •• • . \ .. v.v >* . 

.l-'Cv ‘A*’. ‘v, - *-v’- V.L 

*-* ; ’ : - ' ’ • . •’ • * I t v '• 


TdI7)' 


CclS) ' -pi:xr7: 


" •••• • : .o • ' /•’» /. ;• A 

; Ji&j .do 1". 2" ;v ''Vr-i • ^5: . XIPliSSESK 

x 2 - - - "a 2 ' ' ~ ^ 

.• • ^ .J .V -^v. . , y *0*; V V "■* V 

i • • ‘ !‘‘j; , •.5;*.* ; '. r !** . . J • ' .* l *. ” * + *• v • 


dal • 


«, ,,‘V . . * / 


l 


i2 -4- 1 * *2" f il lll 

c:t dt' 


< c ; •? ,* : v r t ( ci 1£»2> J 


( d : •/ ) 

i. 

f 

( c 2 v > 


ORIGINAL PAGE IS i 2 
OF POOR QUALITY 


( c 21) : c 1 * r- 1 i d C: 2 1 r- 2 ■- 1 i 


i 2 ( > 

r’t 


(p22) :■ r t ;r,d C 2 ) ? " 


dal 2 2 

12 a2 


( d 2 2 ) 

X 




da 1 

2 d c< .1 

dt 

<i2 c 

2 f il ) - 


dt 

d t 

0 

da 2 2 

dal 2 2 

c ci : 

(---> 

i 2 <---) C: 2 

'il • 

dt 

dt 

dt 

O 






dal 

r- 1 


C 

- . 3 


dt 

O 


.. 

12 a 2 + il- 


) n 2 ? 

d 2 2 ) t 


• dr-2 2 

do 1-2“ 2 


> ( ) 

12 ( ) a2 

i 1 < - 


( c 2 i ) 




d t d V 

o '■» 


a iff (h.cl > } 


2(1+12 phi2 2 ) 


60 


I 

I 

I 

I 

1 

1 

I 

I 

4 

■ 

J 

J 

4 

I 


1 

w 

I 

I 

I 


<d 2 ?) 

(c 23) dr-'l-ri; 

( d 2 8 ) 

Cc2?> d i f f < h / os 2 ) r 
( ci 2? ) 

( czo ) Jr?.:;'; 

( d 3 0 ) 

(c 21 ) d 2 :■ * ( -1 ) ; 

< d 3 1 > 

(c 22 ) d 2 l 7 . i 
( d 3 2 :* 


ORIGINAL PAGE IS 
OF POOR QUALITY 


dr-1 - 0 


del 2 

i2 c<2 

dt 


del 2 

i 2 ( ) a 2 

dt 


dal 2 

- i.2 ( ) e2 

dt 


dal 2 

- 12 (— -) , a 
d t 


(c33) dependc<Cg3»Cal»a2rolphol»3lphs23> f • 

Cw(el> a 2 7 olF-hfJl .» 3lPho2)D 

<c34) diff(wral)» 


< d 34 ) 


(c 35 .) F 10 17 .} 


< d35 ) 


(c3o) diff <WfQ2) f 


<d3 6) 


<c3?) f2Q;Y.} 


( d 2 7 ) 


dw 

del 

" d w "**; 
del 

"dw 

tm MM M* 

da2 

dw ' 

dq2 


61 


ri 




( c 2 9 ) alf*ha2JplO''2J 


(d 28 > 


du 2 

( ) 

' dal 


( c 3 v ) alPhal*clPh32/<2*C l+a2''2> >+p20'*2/2+v? 


( d 3 v ) 


du 2 dw 2 

do 2 "dol *■-" 

c 1 p h a 1 ~ + + S' 




2 ( q 2 ' + 1 ) 


( O'! 0 ) £ O 1 V C ( V. t F 2 0 ) J 


dw 2 2 2 

5 <&rt(- < ) - 2 o.2 v - 2 v + -2 olphol a 2 + 2 itlPhaUr • 

du d«i ; v . w r:: 

(d 4 o> •: — = - — t— - — -r— ~~ — r 

UU J - ■ 

sar - tCa2 ' 

dw 2 2 ■; ’ 2 

sort (- — ) - 2 g 2 v - 2 v + 2 slphcl • a2" " ”+ 2 “al'phsl 

dw dal : 


sc<rt(a2 +1) 


( c 4 1 ) p a r t ( d A 0 » 2 ) i 


( d •) 1 ) 


dw 2 2 2 <'•' 

sart <- ( — -> - '2 a 2 v 2 v + 2 dlr-hdl a2 +' 2 alpftoV) 

dw dal ' , 


cart<a2 + 1) 


( d * 1 2 ) 


C c *{ 3 > pcrt(d41»2)? 


dw 2 

cf (aiphcir r. )2 : ; 

dal ‘ \ : 


M'- ' 

. ! rl/" '■ :. s • • - 




dw 2 


<d43> 


sort < - < — :2 a2 •• v - 2 *y. t 2 -alpha 1 ,- a2 r + 2 "a. 1 Ph ol^7" p '• 1 [ W: 


;art<a2 + 1) 


ORIGINAL PAGE IS 
OF POOR. QUALITY 


62 



^ i. 


I 

4 

1 

I 


( c -1 ’ ) ' i nted r 3 1 o ( 2 ; a2 ) 5 


■ORIGINAL PAGE IS • i 
OF POOR QUALITY 




dw 


/ s a r t < - ( o - 2' a ° v - a 


(dAA ) 


L' . dal 

I 


2 a 2 V - 2 V + 2 alphal q2 h 2 alpha! ) 


- - - - - "dd 2' 


4 / sart(&2 + 1) 

f 

(cA?) w*clphs 2 ' % ( l/2)*al+d44F<H2+Kf 


du 2 


/ sort <~ ( ) - 2 al v - 2 v +’2 alphal a 2 + 2 alphdirr 

C dal 


( d *1 7 ) w =-- Cl 


2 

/ 


d«2.' 


sc«rt(a2 + 1) 


! dw ! * '""'dw 

+ al ! ! l; + f< alpha If <---') 

Idol ! dal 


( cM S ) 
< d A 2 > 
( e A ? ) 


s -• 2 # ( a 1 p h s 1 - v ) ! 


3-2 (alphal - v> 


c-alphs2/ ( 2* ( alpha 1-v ) ) ? 





( d v > 


dw 2 

( ) 

dal 


2 (alphal - v) 


<c51 ) ' intcsr'3te< 1/C l-c/< l+a2"2> >"< 1/2) »a2> » 


/ 

C 1 

( d E 1 5 ' I :- : do2*‘ 

2 c 

/ sartCl —)'-•• 

a 2 + 1 







; • . 5 ***‘ ' 



( cS2) dif f ( w » 2 1j- hsl ) = l/3"< l/2)$2» 


/ 

C 


1 ■ 


da 2 


'{d52> 

la S3) .« o rt < d52 r 2 ) ? 


clu 

ds i pita 1 


/ sart< 1 - 

£_ : _, 

sort <3 ) 


/ 

i 

I -• 

3 

/ cert(l - 


da 2 


) . 


( d53) 


a2 + 1 

. ,o o r t ( 3 > 

* ~~ > ' lr,ia - r - t ^< l/(l-c/(l + u2"2) )-<l/2) rQ2»07C«20> J 

a 20 


< r j 2 ) 


C 

X 

3 


' 1 •• . 

dc .,, 

c 

ccrt(l ) 

o 

o2 + i 


' t»cf to i = < 1/3" < 1/2) 


a 2 0 


/ 

c 

I 

3 

/ 


( d 3 ) 


bet si ~ 


scirKl ■ — ) 

- •* ’ • r> ■ 

, • q 2 , + l 


d c< 2 


sort Co J':.": 


' °** GTNal page is 

OE POpR Quality 


64 


111 <■! / (C H-o2"2) * < 1 -c/( l +o2"2) ) "( lVs, ) , ct2 ; 0 ra20 > ; 


a 20 


M-3 ) 


C 

I 

:i 


- lot 2' 


(o2 i i) C C5 P t ( 1 


) 


Cl 2 + 3. 


<eS) bots2=(i/(2tBl.»lian-( 1 /!)).* (olOK clFhal" a/ 2 ) / c '( 


1/2) )#X> 


T ., tS 


ci 20 


/ 

L~ 


s&rUolr-hi‘l) I 
1 


dci2 


/ M2 +1) septa 


a2 + 1 


i 


MS) &fcte2 - 


... + olO 

cart(;) 


nciri<2) c&rt (olf»hsl ) 

(c6) b a t z 2 -- ( 1 l/<2$ si : i h<3l~( 1/2 ) ) ) * ( «10+< slrhsl" <l/2)/a n l 1/2) ) *cM) ? 

a 20 




/ 

r 


eartCalphsl) I 
3 


da 


0 


M2 + 1) 2 a r t ( 1 - ) 

o • 


a 2 i 1 


M6> beta2 - 
<c?) d2; 


sprt < c ) 

2 6Cirt<olr-h'*l) 


+ Q 1 0 


ci 20 

/ 


( d 7 ) 


C 

'I 

3 

/ 


1 


s ci r t ( 1 


--w- ; ;(|q2 

c 

• > 


• r . 


Cl 2 + 1 


(c3) Gi=bets2+3lpha2^*<; 


t j 


( dC ) 


. 1 
(.! 1 ~ o 1 r- h c 2 I 

2 

/ cctrt(l - ■ 
0 


A '■) v <. -. ••> 

to. I u* 1 j U <. 


q2 + 1 


( C ? ) Ci 5 ; 


C ‘ 1 

sort ( blrhsl > I * — 

* . 3. 2 

/ <a2 + i > cart(l 

o 


c a r t ( s ) 


c<2 * 1 — 

+ 


(d?) bets2 -- 


2 ccrt(clfhol) 


( clO ) facitc2-’dw/dsL y h b 2 i 


(dlO) 


b o t o 2 = — - 

d s 1 >• h z 2 


(ell) 'intosirotot (l/2>*< i/<o-b/<l+a2~2> ) > ~ ( 1/2 ) * < ( 2#sl r--h»2 > / a + o2"2 > ) r «2) ; 


(dll) 


c l 

Clph,;2 I ' dci2 

?, 2 b 

/ <c*2 + 1) o&rt<3 > 


( c 1 2 ) - XX. ( - 1 ) ; 


c t 2 + 1 


< d 1 2 > 


/ 

C ’ 1 

- c 1 r- h b 2 I 

2 2 ■ - ' ' b 

7 (c<2 +1)' sartX's' - ) 


ORIGINAL PAGE IS. 
QE .POOR QUALITY 



i 


(cl 3) bet c2 = c 1 F.’i ; 


ORIGINAL' PAGE IS 
OF POOR QUALITY 





/ 

C 

(dl^!) Cel - I 

II 2 

/' cj2 sq rt ( 3 


CclS) closofilsOf 


1 

b 

) + 'sart(a - 

r» 

c<2 FI 



b 

n 

c2 F 1 


da2 + beta23 



References: 


1. Goldstein, H., Classical Mechanics, Addison-Wesley, Reading, 

Massachusetts, 1951. 

2. Singh, S.N. and Schy, A. A., Invertibility and Robust Nonlinear Control 

of Robotic Systems, presented at 23rd IEEE Conference on Decision 
and Control, Las Vegas, Nevada, December 12-14, 1384. 

3. Singh, S.N. and Schy, A. A., Robust Torque Control of an Elastic Robotic 

Arm based on Invertibility and Feedback Stabilization, preasent at 
24th IEEE Conference on Dicision and Control, Fort Lauderdale, 
Florida, December 11-13, 1985. 

4. Punjabi, A. and Divins, J., Kinematics of Robot Arm Using MACSYMA, 

presented at NASA HBCU Seminar, April 20-22, 1986. 

5. Punjabi, A. and Bivins, J., Solutions of Hami lton-Jacobi Equation for 

Robot Arm Using Symbolic Manipulation Program, presented at Virginia 
Academy of Sciences Annual Meeting, May 16, 1986. 


68 



FIGURES: 



Fig. 1. Single link perpendicular robot arm. <t>, is the angle of rotation 
of first arm and <t >2 is the angle of rotation of second arm. 
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t 

ig. 3. Analytical solution of H-J equation with initial conditions 
<t>1 (0) = 0, ^(O) = 0, p-|(0) = 1, p 2 (0 ) = 1, and V = 1. 
y axis represents $2 (in dimensionless form) as a function 
of time, and x axis represents time in seconds. 
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KINEMATICS OF ROBOT ARM USING MACSYMA 


A. Punjabi and J. Bivins 
Department of Mathematics 
Hampton University 
Hampton, V a 23668 

ABSTRACT 


The Hamiltonian of the robot arm shown in the accompanying figure is 
given by: 


H 


P 


2 

1 


2(I 1 +I 2 ^ 2 2) 


+ 



+ V 


We obtain the Hamiltonian and attempt to solve the Hami lton-Jacobi 
equation for the characteristic function F We also attempt to get the 
Poincare surface of section for this problem. All this is done using one of 
the Symbolic Manipulation Programs. Our final goal is the study of 
nonlinear dynamical phenomena arising in this system. 
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SOLUTION OF THE HAMILTON-JACOBI EQUATION FOR SINGLE LINK ROBOT ARM USING 

SYMBOLIC MANIPULATION LANGUAGES 

Alkesh Punjabi and James Bivins 
Department of Mathematics 
Hampton University 
Hampton, VA 23668 

First the Lagrangian L for single link robot arm is set up in the fixed 
system of coordinates. Then the Hamiltonian H and the Hamilton's equation 
of motion are obtained. Hami 1 ton- Jacobi equation for the characteri stic 
function W is set up. H-J equation is solved in closed form. All these 
steps are done using the symbolic manipulation languages MACSYMA and SMP. 
(Supported by NASA under MACSYMA project.) 
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